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ABSTRACT 

We derive the electric and magnetic form factors of the neutron in the framework of a 
relativistic constituent quark model. Our parameter free prediction agrees well with a 
recent, accurate measurement. The relativistic features of the model and the specific 
form of the wave function are essential for the result. Comparisons are made to other 
models based on VMD, PQCD and QCD sum rules. 

A recent measurement 1 of the neutron electromagnetic form factors, Geu{Q 2 ) and 
Gmu(Q 2 ), greatly increased the Q 2 range of previous data 2 and has significantly 
smaller errors. For the first time it is therefore possible to distinguish theoretical 
models with respect to experimental data from the neutron form factors. In the 
low Q 2 region, vector meson dominance (VMD) models 3 are traditionally used 
to make predictions for the form factor. For sufficiently high momentum transfer 
perturbative QCD (PQCD) 4 predicts the Q 2 dependence of the form factors. To 
describe the behavior at intermediate values of Q 2 the parameterization of Rcf. 5 
uses the VMD form at low Q 2 , constrained by PQCD results at high Q 2 . There are 
additional models which predict the neutron form factors. Reference 6 describes a 
relativistic constituent quark model which is similar to our approach. QCD sum 
rules are used in Ref. 7 to fix the parameters of the soft quark functions for calcu- 
lating the form factors. None of these theoretical models are in good agreements 
with the data for both form factors. 

We recently investigated the predictive power of a relativistic constituent quark 
model formulated on the light-front. 8,9 It provides a simple model wherein we have 
overall an excellent and consistent picture of the magnetic moments and the semilep- 
tonic decays of the baryon octet. The parameters of the model have been fixed in 
Ref. 9 so that we have a parameter free prediction of the neutron electromagnetic 
form factors. 

The light-front dynamic is a convenient scheme for dealing with a relativistic 
system. If we introduce the light-front variables p^ = p° ±p 3 , the Einstein mass 
relation p^p^ 1 — m is linear in p~ and linear in p + , in contrast to the quadratic 
form in p° and p in the usual dynamical scheme. A consequence is a single solution 
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of the mass shell relation in terms of p , in contrast to two solutions for p° : 

p- = (pl+m 2 )/p+ , p° = ±^p 2 + m 2 . (1) 

The quadratic relation of p~ a,ndp± = (p 1 ,^ 2 ) in the above Equation resembles the 
nonrelativistic scheme, 10 and the variable p + plays the role of "mass" in this non- 
relativistic analogy. It is therefore a good idea to introduce relative variables like 
the Jacobi momenta when dealing with several particles. As in the nonrelativistic 
scheme such variables allow us to decouple the center of mass motion from the inter- 
nal dynamics. The light-front scheme shows another attractive feature that it has 
in common with the infinite momentum technique. 11 In terms of the old fashioned 
perturbation theory, the diagrams with quarks created out of or annihilated into 
the vacuum do not contribute. The usual qqq quark structure is therefore conserved 
as in the nonrelativistic theory. It is, however, harder to get the hadron states to 
be eigenfunctions of the spin operator. 12 

The light-front formalism is specified by the invariant hypersurface x + = x° + 
x 3 = constant. The following notation is used: The four- vector is given by x — 
(x + , x~ , x±), where x ± = x° ± x 3 and x± = (x 1 ^ 2 ). Light-front vectors are de- 
noted by an arrow x = (x + ,x±), and they are covariant under kinematic Lorentz 
transformations. 13 The three momenta pi of the quarks can be transformed to the to- 
tal and relative momenta to facilitate the separation of the center of mass motion: 14 

P = P1+P2+P3, 

p El n - pt+pt 

k± = (1 - Of>l-L - &2±, 
K± = (1 - V){Pl± +P2±) ~ VP3.L- ( 2 ) 

Note that the four- vectors are not conserved, i.e., P1+P2+P3 7^ P- In the light-front 
dynamics the Hamiltonian takes the form 

where M is the mass operator with the interaction term W 

M = M + IT, 

wa K ± M i ™ 2 

M o = ft \ + + 1 ' 

r?(l-r?) T) l-v 

M| = |±f , (4) 

with m being the mass of the constituent quarks. To get a clearer picture of M we 
transform to fc 3 and K 3 by 

= Ei +k 3 = E 12 + K 3 
Ei + E2 E\2 + E 3 

E 1/2 = (k 2 +m 2 /2 ) 1 / 2 ,i?3 = (K 2 +™ 2 ) 1 / 2 ,i? 12 = (K 2 + M 2 ) 1 / 2 , (5) 



The Electric and Magnetic Form Factors of the neutron 3 



where k = (fci, &2, £3), and K = {K\, K 2 , K3). The expression for the mass operator 
is now simply 

M = £12 + E 3 , M 3 = £1 + £ 2 • (6) 
The diagrammatic approach to light-front theory is well known. 15,16 It provides 
in principal a framework for a systematic treatment of higher-order gluon exchange. 
In this work we limit ourselves to the tree graph. Since we set Q + = we can 
preserve the correct qqq structure of the vertex. All relevant matrix elements that 
we investigate arc related to 



(p \qi + q\ p) Vp'+p+ = m+, (7) 

where the state \p) = \p) / \fp~^ is normalized according to 

(p'\p)=S(p'-p). (8) 
The matrix element M + can be written in terms of wave functions as: 9 

*' + '^S d ' ki ' K {WiSrS' 2 ^- K ' Wk ' K) ' (9) 

where K' ± = K± + r)Qx, and N c being the number of colors. 

The electromagnetic current matrix element for the transition n — * n 7 can 
be written in terms of two form factors taking into account current and parity 
conservation: 



(n',\'p'\J»\n,\ P ) = ux(p') 



ux(j>) (10) 



with momentum transfer Q — p' — p, and the current J M = eq^^q. In order to use 
Eq. (^|) we express the form factors in terms of the + component of the current: 

Fi(Q 2 ) - (n',T |J + M), 
Q±F 2 (Q 2 ) = -2M„(n',T \J + \n,i)- (H) 

For Q 2 = the form factors Fi and F2 are respectively equal to the charge and 
the anomalous magnetic moment in units e and e/M/v- The Sachs form factors are 
defined as G M =F 1 +F 2 and G E =F l - tF 2 with r = Q 2 /4Af2. 

Since the center of mass motion can be separated from the internal motion, the 
wave function ^ is a function of the relative momenta k and K. The product ^ = 
$X0 with $ = flavor, x = spin, and 4> — momentum distribution, is a symmetric 
function. The neutron wave function is given by: 

* = —= (ddux X3 + permutation) cj), (12) 
V3 



with 



xf = ^(m + UT-2TU), 

xf = ^(2UT-m-TU). (13) 
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Since the wave function $ must be an eigenfunction of j 2 (J being the total spin 
of the neutron) and the longitudinal component j'3, the spins f and J, have to be 
rotated by the Melosh transform. 6,12 The S-state orbital function 4>{Mq) is chosen 
to be 

N 

6(Mo) = -1-^ 57- , (14) 

with a and n being phenomenological parameters and N being the normalization 
given by: 

V f d 3 kd 3 K<f> 2 = l. (15) 



The form factors are calculated by inserting Eq. (|l2} ) into Eq. (|J). The result is 
rather lengthy and the explicit expressions are given in Ref. 8. The exponent n is 
fitted to the proton form factor Gup giving n — 3.5. The constituent quark mass 
m and the length scale parameter a are fitted to the proton magnetic moment and 
the weak neutron decay, which results in m = 0.263 GeV and a = 0.607 GeV. 



Fig. 1. The magnetic form factor of the neutron compared with the dipol fit, G Mn / ' fJ-nG n ■ The 
experimental data are taken from Ref. 1 with statistical and systematical errors. Solid line, our 
calculation with pole type wave functions; dashed line, our calculation with a harmonic oscillator 
type wave function; dash-dotted line, VMD model from Hohler; dash-double-dotted line, Gari- 
Kriimpelmann model; dotted line, QCD sum rule prediction by Radyushkin. 



Figures 1 and 2 show the magnetic and electric form factors of the neutron re- 
spectively. The figures give the deviation from the dipol fit Gd = (1 + Q 2 /My)~ 2 
with My = 0.84 GeV. Only experimental data from SLAC NE11 1 are given since 
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Fig. 1. The electric form factor of the neutron compared with the dipol fit, Gi / G 2 -). The data 
and curves are marked the same as in Fig. 1. 



previous data do not distinguish between the various theoretical predictions. The 
present calculation (solid curves) is in very good agreement for both form factors. 
There is only a slight deviation for the magnetic form factor around 2 GeV 2 . To 
show that the specific form of the wave function in Eq. is essential for the 
result we compare the result with the commonly used exponential wave function 
4>(M ) = N exp (— Af 2 /2a 2 ). We also fixed the parameters by fitting other elec- 
troweak nucleon properties, 9 and get m = 0.267 GeV and a = 0.56 GeV. The 
dashed line shows a rapid decrease for Gmu at already 1 GeV 2 , which indicates 
that the exponential wave function is not useful at that energy range. In the non- 
relativistic limit, a/m — > 0, the form factors fall far below the dipol fit for any 
reasonable value of a and m. a The relativistic treatment is therefore important, 
which is a fact also observed for the pion. 17 The VMD model (dash-dotted curves) 
from Holder 3 agrees with the Geu data, but overestimates Gmu- The model from 
Gari and Kriimpelmann 5 (dash-doublc-dotted curves) predicts F\ n = 0. It is there- 
fore in very poor agreement with Geu, and in addition underestimates Gm-h- The 
QCD sum rule predictions from Radyushkin 7 (dotted curves) agrees for Geu and 
underestimates Gmu, approaching Gmu for high Q 2 . The QCD sum rule is not valid 
in the infrared region Q 2 < 1 GeV 2 due to singular power corrections at Q = 0. b 

a For the wave function in Eq. (|l4|) we get Gmyi = Atn(l+2Q 2 /(a 2 +9m 2 )) - 3 5 as the nonrelativistic 
limit for low Q 2 and G M „ = /i„(v / 2Q/3m) 1 / 2 (2Q 2 /(o 2 +9m 2 ))- 3 5 for high Q 2 . 

new method for QCD sum rules in the infrared region < Q 2 < 1 GeV 2 is described in Ref. 18. 
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We conclude that the precise measurement of form factors at intermediate ener- 
gies gives valuable constraints on theoretical models. We showed that a model, that 
is in excellent agreement with the electroweak properties of the baryon octet, gives 
a parameter free prediction of the nucleon form factors, which is in good agree- 
ment with recent experimental data. The relativistic features of the model and the 
specific form of the wave function are essential for the good result. 
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